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Abstract 

We investigate the exact behavior of the energy density of a real massless scalar field inside a cavity with a sin- 
gle moving mirror executing a resonant oscillatory law of motion, satisfying Dirichlet boundary conditions at finite 
temperature. Our results are compared with those found in literature through analytical approximative methods. 
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1. Introduction 

Since the pioneering work of Moore [ 1 ] and the subsequent works published in the 70s (for instance Refs. Hill), 
the problem of particle creation from vacuum due the interaction of a quantum field with moving mirrors (dynamical 
Casimir effect) has been subject of intense theoretical study (for recent reviews see Refs. [5] and references therein). 
On the experimental side, the dynamical Casimir effect has not yet been observed, despite some experimental schemes 
have been proposed based on simulation of a moving mirror by changing the reflectivity of a semiconductor using 
laser beams [6] or, more recently, using a coplanar waveguide terminated by a superconducting quantum interference 
device |7J|. 

In Ref. 01, Moore considered a real massless scalar field in a two-dimensional spacetime inside a cavity with 
one moving boundary imposing Dirichlet boundary conditions on the field. He obtained the field solution in terms 
of the so-called Moore's equation. Exact analytical solutions for particular laws of motion of the boundary [8], and 
also approximative analytical solutions for the Moore's equation were obtained but, so far, there is no general 

technique to find analytical solutions for it. On the other hand, Cole and Schieve developed a geometrical approach 
to solve numerically and exactly Moore's equation 1 1 ill , obtaining the exact behavior of the energy density in a non- 
stationary cavity considering vacuum as the initial field state 111 llll2ll . The dynamical Casimir effect was also studied 



with different approaches from that adopted by Moore, via perturbative methods for a single mirror [13] and for an 



oscillating cavity 01411 . 

The quantum problem of moving mirrors with initial field states different from vacuum was also analyzed Q 
15 , 16i Fnl Il8ll . Specifically, a thermal bath was investigated for the case of a single mirror lfl5ll . and also for an 
oscillating cavity 111 611180 . In a recent paper [ 19], the present authors investigated the behavior of the energy density 
inside a cavity with a moving mirror, for an arbitrary initial field state, obtaining formulas which enabled us to get 
exact numerical results for the quantum radiation force and for the energy density in a non-static cavity for an arbitrary 
initial field state and law of motion for the moving boundary. However, in Ref. 11911 we applied our formulas only to 
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the case of an expanding cavity, with a non-oscillating movement and relativistic velocities. In the present letter, we 
consider the exact approach developed in Ref. Il9tl to study the behavior of the energy density for a non-static cavity 
with a thermal bath and a resonant oscillatory law of motion. We compare our results with those found by Andreata 
and Dodonov in Ref. 1U8I1 where this problem was investigated through analytical approximative methods for small 
oscillations. We show the limitations in the results obtained via this approach as the amplitude of oscillation grows 
to outside the small oscillations regime. We also investigate the energy density outside this regime for the first time. 
Additionally, we show that for small amplitudes of oscillation our results are in excellent agreement with lfl8ll . 

This Letter is organized as follows: in Sec. |2]we obtain the field solution and the exact formulas for the energy 
density taking the initial field state as a thermal bath and considering that the mirrors impose Dirichlet boundary 
condition on the field. In Sec. [3] we study the behavior of the energy density for the case of an oscillatory law of 
motion for the moving boundary. Finally, in Sec. |4]we summarize our results. 



2. Exact formulas for the energy density 

Let us start considering the real massless scalar field satisfying the Klein-Gordon equation (we assume throughout 
this paper h = c = ks = 1): (d 2 - d 2 ^ji//(t,x) = 0, and obeying Dirichlet boundary conditions imposed at the static 
boundary located at x = 0, and also at the moving boundary's position at x = L(t), that is if/(t,0) = if/(t,L(t)) = 0, 
where L(t) is an arbitrary prescribed law for the moving boundary with L(t < 0) = Lq, with Lq being the length of the 
cavity in the static situation. 

The field in the cavity can be obtained by exploiting the conformal invariance of the Klein-Gordon equation Q1Q]- 
The field operator, solution of the wave equation, can be written as: 

oo 

$(t, x) = 2 [a„(A« (f, x) + H.c] , (1) 

n=\ 

where the field modes if/ n (t, x) are given by: 

iff„(t, X) = p [ <p n (v) + (finiu)] , (2) 

with <p n (z) = exp [-innR(z)], u = t — x,v = t + x. The function R satisfies Moore's equation: 

R[t + L(t)]-R[t-L(t)] = 2. (3) 



Considering the Heisenberg picture, we are interested in the averages (...) taken over a thermal state. For this particular 
field state we have {a n a„>) = 6„„i^„(T) and (a n a„>) = (a„a ,) = 0, where (T) = {exp [nn/ (L(jT)] - 1} and T is the 
temperature. 

Taking the expected value of the energy density operator T = (foo(f, x)), where J2l 



f o0, x) = - 



dtj \dxj 



we can split the renormalized energy density T as follows Jl9ll : 

'T" — ^vac + ^"non-vac i 

where 

Tvac = -/(V)-/(M), 
Tnon-vac = -g(v) - g(u), 



with 



g{z) = -\Yjn[R'{z)fun 



(4) 

(5) 

(6) 
(7) 

(8) 



f( , _ J_ I R^d _ 3 

/U; 24tt I /?'(z) 2 



«"(z) 



#(z) 



(9) 



In Eqs. ([8]) and (O the derivatives, denoted by the primes, are taken with respect to the argument of the function R. 
For further analysis, it is useful to write lfl9ll : 



T = -h(v) - h(u), 



(10) 



where h(z) = f(z) + g(z). 

In the static situation t < 0, where both boundaries are at rest, the function R is given by R{z) = zILM- The 
functions / and g, now relabeled, respectively, as and g^ s \ are now given by: 



48Lg ' 



(11) 



(12) 



n=l 



Note that in the static situation 7~ V ac is the Casimir energy density 7~ cas = -7t/(2ALfy. In Ref. IU9I1 it was shown that 
the behavior of the energy density in a cavity is determined by the function h, which obeys: 



h [t + L (f)] = h [t - L (f)] Mi) + S(t). 



where 



i 



1 - V ( t) 



1 +L'(t) 

L"'(t) 



(13) 



(14) 



12tt[1+L' (t)f[l-L' (t)] 
l L" 2 (t)L'(t) 
4tt[1 +L'(f)] 4 [l -L'(t)f 



(15) 



Eq. (TT3T > enables us to obtain recursively the value of h(z), and consequently the energy density T (TlOl l. in terms of its 
static value 



Solving recursively the Eq. (fT3l . as discussed in Ref. [ 19] in details, we can write h(z) in the following manner: 

h(z) = h (s) Mz) + S{z), (16) 

where: 

_ " (z) 

&k) = Y\wm>\, d7) 



S(z) = ^S[?/z)]]~[^[/ i (z)], 



(18) 



7=1 



f i+] + L(f,- + i) 



fi + Uh), 

t- t - L(ti), i = 1,2,3. 



(19) 
(20) 



being n the number of reflections of the null lines on the moving boundary world line during the recursive process. 
We remark that in Eq. ( TTol l. the functions J[ and S depend only on the law of motion of the moving mirror, whereas 
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the dependence on the initial field state is stored in the static value /z (s) . We also observe that, for a generic law of 
motion, J[ and S are different functions, with the following properties lfl9ll : 

ft(z) > V z , &(z < La) = 1, S(z < L ) = 0, (21) 

The energy densities 7" vac and 7"non-vac (Eqs. (O and (0) are now respectively rewritten as: 

Tvac = [^f(u) + Mv)] - S(u) - B(v), (22) 

Tn 011 . vac = [iu» + ^f(v)] . (23) 

From Eqs. (0, (l22b and ( 1231 the exact formula for the total energy density T is now given by: 

T = -h {s) [ft(u) + iK{v)\ - S(u) - S(v). (24) 

In this two-dimensional model the instantaneous force T acting on the moving boundary (disregarding the con- 
tribution of the field outside the cavity) is given by T(t) = T[t,L(t)]. We also define 'Fvac(f) = Tyac [f , L(t)] and 

^^non-vacW — Ton-vac [f> L(f)]. 

With this results in hand, we are ready to investigate the exact behavior of the energy density for a thermal state 
and make comparison with the results found through the analytical approximative approach found in the literature. 



3. Comparing exact and approximate results for the energy density 

From Eqs. ( TT4"1 >, ( fT5l ), ( TPTb and ( TT~8T > we see that the functions J[ and & are different one from another for an 
arbitrary law of motion. Therefore, our first conclusion is that the functions 7~ V ac and 7~ n0 n-vac (given by Eqs. (1221) and 
(1231) ) consequently have, in general, different structures as well. Hereafter we use the word structure in the following 
sense: two graphs have the same structure if they have the same number of maximum and minimum points and these 
points are at the same positions in both graphs. As a direct consequence of our first conclusion, we can say that the 
thermal force 'Fnon-vac and the radiation reaction force T\ ac have in general different structures. At a first glance, our 
conclusion contrasts with that found by Andreata and Dodonov lfl8ll . according to which T vac and 7" n0 n-vac exhibit a 
same structure for initial diagonal states (as the case of the thermal state). Next we will discuss this issue. 

Let us consider the particular laws of motion given by 



L(t) = Lo 



1 + e sin 



put 



(25) 



where Lo = 1, P = 1.2, and e is a dimensionless parameter. This oscillatory boundary motion was investigated in 
several papers (see, for instance, Refs. Jgl [T§]), with the calculation of the energy density developed in the context of 
analytical approximative methods, considering small amplitudes of oscillation (|e| <K 1). Taking as basis the results 
found in Ref. 11811 . the renormalized energy density T, corresponding to the laws of motion (f25j) is given by T « 7~ (a ', 
with 

7-(«) = _ ( ftM _ p 2 f^)[ s ( u ) + S ( V )] - 2//«, (26) 



where: 

s(z) 



.2x2 



[1 + k 2 + 2(-1)Pkcos(pttz)] 2, 
sinh(/?r) 



1 + sinh 2 (pr) 



1 

t - ent, 

2L Q 



being t - Nip and a non-negative integer. The energy density T (a) in Eq. d26l i also can be written in terms of the 

r< fl ) 

non-vac > 

Tvac = (P 2 ~ D/ (S) W") + *(v)] - 2p 2 f s \ (27) 



vacuum and non-vacuum contributions as T (fl) = T^ffl + < Tj5L.- ms ., where: 
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Tnon-vac = -£ ( ' 5) M") + sfy)l (28) 
Let us focus initially on the case p > 1. Since (p 2 - l)/ (s) > and -g (s) > 0, Eqs. (f27T > and (1281 give that the functions 
7~iac an d ^non-vac have the same structure 1 1 8] . To conciliate this conclusion with our first conclusion mentioned 



above, we will show next, starting from our exact approach, that we can find a class of motions for which the energy 
density (l24l i is T ~ 7~ (fl) , so that TVac and 7" n0 n-vac exhibit approximately the same structure, and that the particular 
laws of motion d25l - investigated in Ref. 11811 - belong to this class. 

Looking for conditions under which T V ac and 7" n0 n-vac have the same structure we find that one condition is 
provided by the laws of motion for which 3\.(z) and S(z) have a linear relation of the form: 

S(z) = krM{z) + k 2 , (29) 

where k\ and k 2 are constants. From the properties given in Eq. (f2TJ, we get k\ - —hi, resulting in: 

fl(z) = k x [&(z) - l] , (30) 

and from Eqs. ( 1221 . ( 1231 and ( l30l . we have: 

Tvac = - (/ <s) + h) [Mu) + My)] + 2*i, (31) 

Tnon-vac = -g W [#(«) + #(v)] , (32) 

T-- (h (s) + *i) \Mu) + #(v)] + 2*i. (33) 

Then, our second conclusion is that if the constant factors multiplying J[{u) + Jl(y) are different from zero and have 
the same sign, as the ratio 

ota-^*- (34) 
Mz) - 1 

becomes more close to a constant value k\, that means 

o-(z)**i, (35) 

more the structures of the functions T V ac and 7" n0 n-vac become similar to one another, whereas if they have different 
sign then where we find valleys and peaks in a graph, we can have respectively peaks and valleys in the other. As a 
direct consequence, for the class of motions obeying the condition (135b the quantum radiation forces iTvac and !Fnon-vac 
also have similar structures. 

Comparing our Eq. (l33l with ( 1261 . we see that both have the same structure. Our third conclusion is that the 
formula ( 1261 found by Andreata and Dodonov belongs to the particular class of formulas (given by Eq. (l33b ) for the 
energy density. We will show that this occurs because for the laws of motion (l25l . as the condition |e] «c 1 is better 
satisfied, better is satisfied the condition ( l35l . 

In Fig. 1, using ([PTT i and ( TT8l . we plot the ratio <x (Eq. d34l l) for p — 2, taking into account three values of 
amplitudes of oscillation: e = 1CT 3 , e = 10~ 2 and e = 10 l . We observe that <x is more approximately the constant 
value -4/ (5) for e = 1CT 3 (dash-dotted line) than for e = 10 _I (solid line). 

In Fig. 2 we plot <x for the case p — 3, observing that as e becomes smaller <x tends to the value -9/ (?) . Then we 
see that, in the limit e — > we get: 

cr -» h = -P 2 f s) . (36) 

In Fig. 3, we plot the ratio A(z) = M(z)/s(z) for p — 2, taking into account the values e = 1CT 4 , e - 10~ 3 , e = 10~ 2 
and e = 10~'. We observe that as e — > we have A(z) — > 1, what means: 

Wti -» s(z). (37) 
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-0.5 1 1 1 1 1 1 1 

1 2 3 4 5 6 7 

Z 

Figure 1: The ratio cr{z) = fife)/ ~ *] 

for the law of motion given by Eq. (f25T > with p = 2. We use different 
scales for cr(z) in each case. The solid line corresponds to the case e = 10 -1 . The dashed line corresponds to the 
case e = 1CT 2 , exhibiting 40 x [cr(z) + 4/ (s) ] - 4/ (!) . The dash-dotted line corresponds to the case e = 1CT 3 , showing 
200 x [o-(z) + 4/ w ] - 4/ (s) . The dotted line corresponds to the case e = 0. As e -> we have cr(z) -> -4/ (s) . 




Figure 2: The ratio cr(z) = S(z)/ [j?l(z) - lj for the law of motion given by Eq. ( l25T l with p = 3. We use different 
scales for cr(z) in each case. The solid line corresponds to the case e = 10 -1 . The dashed line corresponds to the 
case e = 10~ 2 , exhibiting 40 x [cr(z) + 9/ <s) ] - 9/ (!) . The dash-dotted line corresponds to the case e = 10 -3 , showing 
300 x [o-(z) + 9/ w ] - 9/ (s) . The dotted line corresponds to the case e = 0. As e -> we have cr(z) -> -9/ ( ' 5) . 



Eqs. ( f36b and d37b complete the mapping between Eqs. ( 1331 ) and d26b . demonstrating a perfect agreement between 
two completely different approaches to the problem 

Now, let us investigate the following point: since Eqs. (l27l > and (l28l i are approximations, we should find differences 
between the structures of the vacuum and non-vacuum parts when we consider the case of the law of motion d25l with 
our exact approach. To investigate this issue, we study the law of motion ( |25] | for e = 10~ 2 and e = 10 _1 . Our aim 
now is to verify, using the exact approach, the similarities and differences between the structures of T V ac and 7" non -vac 
for both values of e. 

Although the formulas ( fTTT i. ( TT~8T >. (l22l and ( 1231 are formally exact, to extract numerical values for 7~ V ac and 
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Figure 3: The ratio A(z) = Jl(z)/s(z) for the law of motion given by Eq. d25l ) with p - 2. We use different scales 
for in each case. The solid line describes the case e = 10 l , exhibited as 1/6 x [A(z) - 1] + 1. The dashed line 
describes the case e = 10~ 2 . The dash-dotted line corresponds to the case e = 1GT 3 , exhibiting 1 /5 x A(z), whereas the 
dotted line corresponds to the case e = 1CT 4 , showing 20 x A(z). As e — > we have A(z) — > 1. 



'/"non-vac we need to calculate the functions f,(z), which are given by Eqs. $1% and (l20b . For the law of motion d25l > 
(and in general) we can get only numerical solutions of $1% and d20b . For a given value of z, the Eq. ( fT9b can 
be solved numerically and the result for t\ naturally has a certain limited accuracy. When t\ is used in Eq. ( |20T i to 
calculate t2, the solution of ?2 + £fe) = h — L(fi) can give a result less accurate than the result previously obtained 
for ?i, and successive calculations of the remaining values of f, via equations ( l20b (for z = 2,3,4...) could generate 
a final result for t n with a poor accuracy, if compared to the accuracy of the initial value for t \ . Moreover, when we 
insert the numerical values for f, in Eqs. (1171 1. (1181 1. (1221) and ( f23l , the final values for 7~ V ac and T n0 n-vac could have 
their accuracy diminished even more. To deal with this question and control the final accuracy of our results, we 
perform the numerical calculations in Maple computer algebra system 12011 . which enables us to control the number of 
digits used when calculating with floating-point numbers. To obtain a final numerical value for TVac or 7~ n0 n-vac, we 
carry out several independent calculations using our routines developed in Maple 112 111 . In each calculation we take all 
numerical solutions with a certain number of digits. Considering from 3 to 100 digits, we observe the convergence of 
the results as the number of digits (related to the initial accuracy considered for the solution of ( fT9b ) is enhanced. This 
enable us to point which are the exact digits (the accuracy) in our results. For instance, the values of 7~ V ac(10, 0.5) 
(L () = l tP = 2) performed with 3, 4, 5, 6, 10 and 20 digits are given, respectively by: 0.856, 0.8638, 0.86364, 
0.863704, 0.8637005768 and 0.86370057587773139184. Analyzing also the convergence of the results up to 100 
digits, we can obtain a final accuracy of 10 or more digits, but we just display the result as 7~ V ac(10, 0.5) ~ 0.864, 
where the first two digits can be considered as exact digits. Hereafter, the exhibited results have accuracy at least up 
to the penultimate digit shown. 

In Fig. 4, using the formulas d22l and (1231 we plot the energy densities T V ac and 7~ n0 n-vac for the case p = 2, T = 1, 
e = 10~ 2 at t — 20.2. We see that both energy densities have two peaks (in this case, located at x — 0.30 and x = 0.70, 
with values r vac (20.2, 0.30) = r vac (20.2, 0.70) ~ 1.94 and r non -vac(20.2, 0.30) = r nO n-vac(20.2, 0.70)~ 23.3), both 
have three minimum points (located at x = 0, x — 0.50, x — 1.01, with values 7~ V ac(20.2, 0)~ -0.453,7~ V ac(20.2, 0.5) * 
-0.322, r vac (20.2, 1.01)* -0.461, r non . va c(20.2, 0)* 0.647, r non -vac(20.2, 0.5) * 1.9, r non . va c(20.2, 1.01) * 0.579), 
so that Tvac and 7~ n0 n-vac exhibit the same structure, as predicted by the approximate analytical formulas ( l27l ) and d28l > 
1 1811 . which are based on the assumption |e| <K 1. We remark that x * 1.01 corresponds to the position of the right 
mirror when t = 20.2 (L(20.2) ~ 1.01). 

For the case e = 10~\ we get for 7~vac(20.2, x) and 7" non - vac (20.2, x) that both present again two narrow peaks, 
located at x = 0.30 and x = 0.70 (similarly to the case showed in Fig. 4, but with values 7~vac(20.2, 0.30) = 
r vac (20.2,0.70) * 0.249 x 10 13 , r non -vac (20.2, 0.30) = r non -vac (20.2, 0.70) ~ 0.126 x 10 14 ). Moreover, in the 
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Figure 4: The energy densities 7~ vac (dashed line) and 7~ n0 n-vac (solid line) at the instant t = 20.2 with T = 1, plotted 
via the exact formulas d22l and 0231 . We consider the law of motion d25T > with p — 2 and e = 10" 2 . 



case e — 10 , 7~vac(20.2, x) exhibits several other maximum and minimum points that are not visualized in the graph 

Of 7~non-vac(20.2,x). 

I 1 1 1 1 1 1 1 1 1 




Figure 5: Detail of the energy density 7~ V ac at the instant t - 20.2, for the law of motion d25T l with p = 2, T = 1 and 
e = 10 _1 , plotted via the exact formula (l22l in the region 0.2 < x < 0.3. In detail we show 7~ V ac in the sub-region 
0.26 < x < 0.28. The spacing used between the calculated points in the graph is 10~ 4 . 



In Fig. 5 we investigate details of the behavior of T vac for the case e — 10 1 in the region 0.2 < x < 0.3. We can 
see peaks for the following values: x * 0.250 and x * 0.289; we also see minimum values at the points x ~ 0.230, 
x a 0.28 and x as 0.295. 

In Fig. 6 we see the behavior of 7~ n0 n-vac (for e = 10~' and T = 1), in the same region as Fig. 5, but there is 
no peak or valley. Then we verify that when we consider e = 10 -2 and e = 10 _1 , since the the former value is in 
better agreement with the conditions |e] <K 1 and cr as -p 2 / (s) = -An, no difference between the structures of 7" V ac 
and Tnon -vac is perceived, but for the latter value of 6 differences between the structures of *7~vac (Fig. 6) and 7~non-vac 
(Fig. 5) become evident, as predicted via the exact formulas (l22l and (l23l . 

For p = 1, Eq. d27l > gives 7~ vac as Ty^. = -2/ (s) = T C as, so that the energy density would conserve its vacuum 
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X 

Figure 6: The energy density Tnon-vac at the instant t = 20.2, for the law of motion d25l > with p = 2 and e = 10" 1 , 
plotted via the exact formula ( 1231 in the region 0.2 < x < 0.3. In detail we show Tnon-vac in the sub-region 0.26 < x < 
0.28. The spacing used between the calculated points in the graph is 10~ 4 . 
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Figure 7: The time evolution of the energy density T va c for the law of motion (l25l . with p — 2 and e = 10~ 2 , plotted 
via the exact formula d22l . The dashed line shows two peaks at t — 19.9, initializing the merging process. The solid 
line shows the maximum peak formed at x — 0.5 and t = 20. The dotted line shows the energy density at t — 20.2, 
after the merging. 



(Casimir) value, whereas Eq. ( l28l gives for 7"n 0n _ vac a spacetime dependence. However, the exact behavior is given 
by Eqs. (l22l and (1231 . which show that the values of both T va c and Tnon-vac change in time and space. For instance, 
for the law of motion d25l l with p — 1 and e = 10~ 2 , the exact behavior of 7^(100.5, x) exhibits two minimum points 
surrounding a peak located at x — 0.5 (see also Ref. lITlll ). whereas 7^on-vac( 100.5, x) displays a peak at x — 0.5. For 
this case p = 1, Eq. (l36l l remains valid, so that as e — > 0, Ty ac and T^ n . vac display more similar structures. 

Now, we will compare exact and approximate results in the prediction of the maximum value of the peaks in the 
energy density moving in an oscillating cavity. In this context, let us analyze again the behavior of T vac (t, x) and 
Tnon-vacO, x) for the case p = 2. The two peaks showed in Fig. 4 (for T vac (t, x) or 7~ n0 n-vac(f, ■*)) have the same value 
(see values presented above), move in opposite direction and at (f, x) = (N, 0.5) (N is a non-negative integer) they 
merge forming a single maximum peak (see Fig. 7). The maximum value of the energy density occurs when the 
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two peaks merge. This value is represented for vacuum and non-vacuum parts, respectively, by T™* and 7^? vac . 
Exactly, we have 

T™ x =T mc (N,0.5), (38) 



* in 



T n , 



=(#,0.5). 



From an approximate analysis, taking as basis the results found in Ref. Ill 8H . T\ 

T (fl)max where 
1 non-vac wnere 



tn f7-max ^ (T-(i)max j ^-max 
Lu ' vac ~ * vac a " u ' non-vac 



max „ n( i rT-max 
vac dllu ' vac 



(39) 

grow in time according 



^-(a)max _ f^ f (s),J,T 
' vac 



6/W( e 8T - l) - 2f 



(a)max 



-2g' 



s) (l+«) 2 
(l-*) 2 ' 



(40) 



(41) 



In Fig. 8 we examine the case given by the law of motion $25[ , p = 2 and e = 10~ 2 , and visualize agreement 
between the growing in time of the peaks predicted by the approximate formula ( l40l i and the exact values obtained 
via (l38l . In Table 1, again considering e = 10~ 2 , we compare T™ x with T^" 3 ^ for larger times, and also compare 
We remark the increasing error of the approximate formulas in comparison with the exact 



T" 

* in 



with T, 



(a)max 



formulas, observing that the error of T^^ c grows more rapidly. 




Figure 8: The energy density T™ x (solid line), given by Eq. (1381 . and its growing in time predicted by (1401 . for the 
case of e — 10~ 2 . 
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Table 1: 7~™ x and T"n™ x vac (for T - 1) computed via numerical exact method (column 2) and via approxi- 
mate analytical formula (column 3), with Ln = 1, e = 10~ 2 and p — 2. The percent error (7~ max , - 

J • ' 1 r x Lvac, non-vacj 

r ( fl )max )/7-max x lQQ j showed in co l umn 4. 

[vac, non-vacj ' ' [vac, non-vacj 



t 


/T-max 
' vac 


^-(a)max 
' vac 


Percent Error 


10 


0.864 


0.856 


0.870 


10 2 


115 


113 


2.16 


5x 10 2 


0.832 X 10 27 


0.761 X 10 27 


8.44 


10 3 


0.175 X 10 55 


0.148 x 10 55 


15.7 




7-max 
' non-vac 


^-(a)max 
' non-vac 




10 


13.12 


13.1 


0.17 


10 2 


109 X 10 4 


107 X 10 4 


1.64 


5x 10 2 


0.785 X 10 28 


0.723 X 10 28 


7.95 


10 3 


0.165 x 10 56 


0.149 X 10 42 


99.9 
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Figure 9: The energy density 7~™ x (solid line), given by Eq. ( l38l l. and its growing in time predicted by (l40l i. for the 
case of e = KT 1 . 



In Fig. 9 we examine the case discussed in Fig. 8, but with a larger amplitude: e = 10~' . Now, we see disagreement 
between the maximum value of the energy density (l38l and its growing in time predicted by (l40b . In Table 2, we 
examine the case discussed in Fig. 9, but for larger instants. We see large discrepancy between approximate and exact 
formulas. 

In counterpart, if we diminish the amplitude to e = 1CT 8 , we get an excellent agreement between approximate and 
exact formulas for short and long times. For instance: 7"™ x (4x 10 4 ) ~ 7"^ max (4x 10 4 ) ~ -0.129; T™ n x vac (4x 10 4 ) ~ 
T" X c(4xl0 4 )*3.75. 
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Table 2: T™ x and T^^ VSIC (for T — 1) computed via numerical exact method (column 2) and via approxi- 
mate analytical formula (column 3), with Ln = 1, e = 10 _1 and p = 2. The percent error (7~ max , - 

J • ' r r Lvac, non-vacj 

r ( fl )max )/7-max x lQQ j showed in co l umn 4. 

[vac, non-vacj ' ' [vac, non-vacj 



( 


rr-max 
' vac 


~-(a)max 
' vac 


Percent EiTor 


10 


191 x 10 4 


113 X 10 3 


94.1 


40 


0.338 x 10 26 


0.266 x 10 22 


99.9 




7-max 
' non-vac 


^-(fl)max 
' non-vac 




10 


971 x 10 4 


107 X 10 4 


88.9 


40 


0.171 x 10 27 


0.252 X 10 23 


99.9 



4. Conclusions 

Considering a thermal bath as the initial field state in a cavity with a moving mirror in a two-dimensional space- 
time, we found for small amplitude of oscillation good agreement between the exact and approximate results for the 
maximum values of the energy densities. This agreement strengthens the validity of the analytical approximate results 
obtained in Ref. Ill 811 . and also reinforce the validity of the exact formulas and numerical results discussed here. How- 
ever, for larger values of e, as shown in Table|2] significant discrepancies appear. This is expected, since the analytical 
formulas are valid for |e| <K 1. Then, we see that the exact formulas (l22l) and d23l can give results for cases of large 
amplitudes, which are out of reach of the perturbative approaches found in the literature. 

We showed that the energy densities 7~ V ac and 7~ n on-vac have, in general, different structures. However, we found 
that these energy densities can exhibit approximately the same structure for a class of laws of motion for which the 
ratio given in Eq. d34l l is approximately a constant value. We also showed that this condition is just satisfied by the 
oscillating laws of motion with small amplitude investigated in the literature, specifically in Ref. 11811 . where the same 
structure for these energy densities was predicted via approximate methods. We verified that for this class of laws 
of motion there is a direct mapping between the approximate analytical formulas for the energy density found in the 
literature and the exact formulas discussed here. On the other hand, we found that for larger amplitudes of oscillation 
the ratio in Eq. d34b becomes far from a constant value, displaying larger oscillatory behavior. This means that that 
7~vac and 7~non-vac can display different structures. Moreover, the exact formulas (l22l i and d23l can say precisely how 
these structures are. Finally, we remark that, beyond the thermal case, the conclusions found in our letter are directly 
extensible to any other initial state whose density matrix is diagonal in the Fock basis. 
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